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ABSTRACT 

The  propagation  of  sound  in  a  laboratory  model 
consisting  of  a  surface  layer  of  water  overlying  a  thick 
bottom  layer  of  sand  was  experimentally  investigated.  The 
variations  of  pressure  amplitude  as  a  function  of  receiver 
depth,  range,  and  transverse  location  were  measured  for  the 
source  at  constant  depth  in  the  surface  layer  of  water. 
Three  frequencies  for  which  only  the  lowest  mode  propagates 
were  used  in  this  investigation.  The  measurements  of  pres- 
sure amplitude  as  a  function  of  depth  and  range  were  found 
in  good  agreement  with  normal  mode  theory,  when  absorption 
in  the  bottom  was  taken  into  consideration. 
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I .  PROPAGATION  OF  SOUND  IN  THE  HOMOGENOUS  LAYERED  MODEL 

The  transmission  of  sound  in  the  ocean  is  dependent  on 
several  environmental  factors.  The  most  important  of  these 
include  the  depth,  the  nature  of  the  bottom,  sound  velocity 
structure  within  the  ocean,  and  the  shape  of  ocean  bottom 
and  surface.  The  purpose  of  this  experiment  was  to  investi- 
gate the  propagation  of  sound  in  a  two  layered  laboratory 
model,  composed  of  sand  and  water.  In  this  study  the  bottom 
and  surface  of  the  model  were  made  flat  and  experimental 
results  were  compared  with  a  simple  theoretical  model. 
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II .  PROPAGATION  OF  SOUND  IN  SHALLOW  WATER 

A.  INTRODUCTION 

The  propagation  of  sound  in  shallow  water  is  a  compli- 
cated problem,  involving  multiple  reflections  of  the  sound 
between  the  surface  and  the  bottom,  so  that  there  is  consid- 
erable phase  interference  among  the  various  multipaths.  At 
relatively  long  ranges,  the  sound  field  is  mainly  due  to  the 
multiple  reflections.  The  transmission  of  sound  in  the 
shallow  water  also  depends  on  several  environmental  factors: 
nature  and  slope  of  the  bottom,  depth  of  water  relative  to 
wave  length,  roughness  of  surface  and  bottom  and  temperature 
gradient.  As  a  consequence  of  the  multiplicity  of  reflec- 
tions from  the  surface  and  the  bottom  and  the  environmental 
factors,  the  predictions  of  sound'  transmission  in  shallow 
water  can  be  much  more  complicated  than  for  deep  water.  Ray 
theory  can  be  used  for  the  sound  field  at  short  range  when 
the  depth  of  the  water  is  relatively  large.  But,  ray  theory 
has  limited  applications.  It  is  inapplicable  in  a  shadow 
zone.  Ray  theory  also  cannot  be  used  for  low  frequencies 
when  the  wavelength  of  sound  becomes  comparable  with  the 
vertical  scale  of  sound  velocity  variation.  Therefore,  in 
practice  one  ultimately  has  to  resort  to  a  wave  solution  of 
the  problem. 

B.  NORMAL  MODES  IN  TWO  HOMOGENOUS  LAYERS 

Normal  mode  theory  provides  an  exact  solution  to  the 
wave  equation  and  is  also  useful  for  dealing  with  depth 
dependent  sound  profile  particularly  those  leading  to 
the  trapping  of  sound.  There  are  various  approaches  to 
the    normal    mode     theory    but    we     will    use    the 
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eigenfunction-expansion  technique,  which  is  relatively 
simple.  In  this  presentation  we  will  follow  mainly  [Ref.  1] 
and  [Ref.  2 ] . 

Consider   the   case   of  two   homogeneous   layers   with  a 
pressure  release  surface  above  as  shown  in  Figure  2.1  . 


Pressure    release 


do 


H 


Z  (z) 


Z   (z) 

in 


c<  .  r, 


c<  .  r. 


Fluid  i 


Fluid  2 


Figure  2.1    Two  homogeneous  layers  with  pressure  release 
surface . 


Let  P  ,C1  and  P ,C    be  the  constant  densities  and  sound 

speeds  in   layer  1   and  2  respectively.  Assume  there   is  a 

point  source  at  z   =  d  in  layer  1.   The  wave  equation  for  a 
homogeneous  medium  can  be  written  as  ■ 
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V  P  -   i   d2B     -    -    a 


where  (Tis  a  source  function  and  is  defined  by 

at 

Where  Q   is  source  strength 

In  our  solution,  we  assume  that  a  point  source  radiating 
a  single  frequency  can  be  described  mathematically  by 


CT  =  4TTAS(r-r)e; 


I  CO  i 


•s 


(2.2) 


Where  5 ( r  -  r.  )  is  the  Dirac  delta  function,  r  is  the 
position  vector.  A  is  the  amplitude  factor  which  is 
assumed  to  be  unity,  and  co  is  the  angular  frequency  of 
the  source.   Solution  to  the  wave  equation  is  of  the  form 


ju>t 


P(r,t)  =   P(r)eJ 

(2.3) 


Substituting  equation   (2.2)   and   equation  (2.3)    into 
equation  (2.1)  we  get  the  helmholtz  equation 


VP  +  kP   =   -4lT6(r-r  )      l2   4, 
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where  k  is  wave  number  defined  by 


k   =   co 
c 


(2.5) 


The  delta  function  for  the  point  source  at  r  =  d  is 


Z(Z-L) 


1   S(r)6(z  -d) 
2irr 


(2.6) 


In  cylindrical  co-ordinates  equation  (2.4)  can  be  written  as 


ri<3  (rd_\,d2        k2   i-P  -   2  6(r)<S(z-d)  ey"' 


(2.7) 


This   is   an   inhomogeneous  equation.     First,    let  us 
consider  the  homogeneous  equation, 


±  d_(rd 


d 

a? 


k   1-P  =     0 


(2.8) 


and  use  the  separation  of  variables  technique. 

We  assume  that  the  sound  energy  is  trapped  in  the  upper 
layer  and  solution  to  equation  (2.8)  can  be  represented  by 
a  set  eigenfunctions . 


P(r,z,t)  =  ey~T 


R(r)Z(z) 


(2.9) 
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Substitute  equation  (2.9)  into  equation  (2.8)  to  get 


[  «     "r  ar^     dr^nndP'rn  n      J 


0 


(2.10) 


or 


L   R     r    d  r  ^    d  r  ' 


+ 


2  2 

r_Li&  +  k 


=     0 


(2.11) 


This  is  only  possible  if  each  part  in  the  square 
brackets  in  equation  (2.11)  is  equal  to  a  constant  hence  it 
can  be  written  as 


R,~dP>  if) 


(2.12) 


(2.13) 


where  kn  is  separation  constant. 

The  boundary  conditions  on  the  depth  functions  are 


ZJO)     =  0 


(2.14) 
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Z  (H)      =  ZJH) 

in  *n 


(2.15) 


UdZ^A 


(2.16) 


lim  Z2„(z)  =  0 


2-CD 


(2.17) 


Rewrite  equation  (2.13)  in  the  form 


or 


dZn 


+ 


&Z.  =  0 


(k  _k„)Z„ 


AZn   +   k,  Z„   =   0 

d  z2 


2/7  -7 


(2.18) 


where 


k   -  /of   if 


(2.19) 


In  order  to  satisfy  the  boundary  condition  the  pressure 
goes  to  zero  when  Z  — >■  oo  ,  the  solution  Z£n  must  decay 
exponentially  with  depth.   The  only  way  that  Zm  can  satisfy 
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the  boundary   conditions   at   the  surface  at   Z  =    H  is   to 
be  sinusoidal  function  of  depth.  Therefore 


Zm     =  AnSm   Kn1  O^Z^H       (2.20) 


Z         -  Be    "  H-E  Z<  oo       (2.21) 


p.  =  k"  '(*) 


<V 


(2.22) 


Waveguide  propagation  requires  that  kn    and  6   must  be  real 
hence 


£±   ^  k   <  Q 
C<    "     C,  (2.23) 


Application  of  boundary  Conditions  at  Z  =  H  results  in 


Ansin  k-  h      -        B  e  (2.24) 

zn         ~  n 


A,ik   cosk-w      =      _i6Be^ 

7)    zn  zn  "q"  trt  n  (2  .25) 
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which   gives 


Bn     -  A„e   -sin  k-H 


(2.26) 


The    depth    function,    therefore,    can   now   be   written   as 


Z.Cz) 


A  sin  k  -z 

n  zn 


0  <  Z    <    fi     (2.27) 


Z(2)      = 


/J  2/7 


H   <    Z 


(2.28) 


These  depth  functions  have  to  be  orthogonal.  Since  there 
is  a  discontinuity  in  the  slope  at  the  boundary  Z  =  H,  we 
need  separate  orthogonality  constant  for  each  layer.  We 
assume  that  the  constant  is  unity  for  the  layer  1  and  *S  for 
the  layer  2 . 


z.(z)   = 


V 


0<  2-sH 


7Z2„(Z)     H<Z<oo 


(2.29) 


In  order  to  evaluate  "2$  multiply  equation  (2.13)  by  Z. 


2       2 


Zmd_4   (k  -k  )Z  Z   = 

w  a  -.2     T       n       n    m 


dz 


0 


(2.30) 
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interchange  the   subscript  m  and   n  in  equation   (2.30)   and 
take  the  difference  of  the  two  equations  to  get 


d_(Z  dZ„   z  dZ^ 


.2      .2 


dz 


=  (kn-km)ZmZn  (2.31) 


Integrate  both  sides  of  equation  (2.31)  from  0  to  °0 


fu2^-1^1--^-^ 


-m  ndZ  (2.32) 


The  orthogonality  condition  requires  that,     for  m  £   n, 
the  R.H.S  of  the  equation  (2.32)  must  be  zero.  Therefore 


d    (Z     dZ,n 
dT^   1*d"z- 


Z1„dZi^dz+T/d_.f2lBidZLW_Z2fldZ.^dz  =  0 

d  2    J  I  dzv      dz  d  z   J  ,„  _ 


(2.33) 


Z1w,dZ,n  _Z„dZ«l  +  7rZ,„dZw     Z1H.dZ„ 
dz  dzJ0         Ldz  dz 


=     0 


(2.34) 


The  application  of  limits  as  well  as  boundary  conditions 
gives 


v-'m^K&i-wK^j  -  o 


(2  .35) 
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This    is    only   possible    if 


_7A     =    o 
A 


or 


7 


A 

A 


(2.36) 


Therefore,  equation  (2.29)  becomes 


f  AflZ(z) 


£  AnZJZ) 


0  <  2  <  H 


H<    Z  <  oo 


(2.37) 


We  will  normalize   the  depth  eigenfunctions  in   order  to 
evaluate  A„,  ,  consider  the  m  =  n  case 


oO 


Z  (z)dz 

n 


(2.38) 


.H2 

A  Z    Z    dz 

n     in     in         + 


A.Z.  Z    dz 


a. 


n     in     2/7 


=        1  (2.39) 


A, 


fVdz  ,    A  fVdz 

0  ^  JH 


(2.40) 
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Now 


.H    2 

Z.  dz 


m 


0 


=       /  i(  1  -  cos  2k -z  )dz 
J  2         zn 


(2.41) 


Integrating  and  substituting  the  limit  we  get 


/ 


H-2 

Z   dz 


in 


=      IT  H L  Sin  k-  H  cos  k    H  1 

2L  kzn  zn  zn     J 


(2.42) 


and 


.OB     2 

Z   dZ 


*o 


W 


■  z1 

H 


(Sin   k7nH-ep^    )e  p"  dz 


(2.43) 


or 


I 


co    2 

2.  dz 


in 


_l_Sin2k   H 
2£  *" 


(2.44) 


Substituting   equation   (2.42)    and   equation  (2.44)    int< 
equation  (2.40)  gives 


2k 


zn 


k   H  -coS(k-H)sin(k -H)  -A&fsinfk   H)tan(k   H) 

zn  zn  in  \Qj  zn  zn 


A 


(2.45) 
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The  depth  eigenfunctions  now  form  a  complete  set  of 
orthonormal  functions. 

Mow  we  have  to  evaluate  the  Green's  function  for  the 
inhomogeneous  wave  equation 


»  "TSrl  drn)+  «    na^+"   ~""UJ    r 


(2.46) 


?^[ldvddf>fe]-?e'.[^(k2-k")z»]=-fs(r,8(z"d) 


(2.47) 


From  equation   (2.13),   the   second  term   in  the   square 
bracket  is  zero;  therefore,  equation  (2.47)  can  be  rewritten 


as 


_  2_S(r)6(z-d) 

r  (2.48) 


Multiply   both   sides   of   equation    (2.48)    by  Zwand 
integrate  from  0  to  co  gives 


y[il_(rdB..)+k2Rirz^dz    __ 


2_S(r)     <S(z-d)Z„dz 

o  (2.49) 


For   m   =    n 


ilfrd&A.   k2,R     - 


2.5(r)Z(d) 
r  " 


(2.50) 
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Solution  to  this  homogeneous  equation  is 


(21 


R(r)  =  -j7TZ(d)H  (k  r) 

~-n  ■>         n  0      n 


(2.51) 


where  for  k-r  >>  1 


H  (  k  r  ) 


■j-     -J<*S-Z) 


1Tk  r 


C2) 


The  Hankel  function  H0  (k^-  r)  describes  the  solution  in 
the  r  direction. 

Therefore,  the  solution  for  the  sound  trapped  in  layer  1 
(0  <  Z  <  H)  is 


PC.z.t)  =  -j 


-,Z 


27T  A'sin(kd)sin(k  z)e;V"'-*"r+T; 

l  r        n  zn  zn 

knr  (2.52) 


Equation  (2.52)  is  invariant   under    the  interchange 

of  Z  and   d   and  is  valid   within   the   region  0  <  Z   <  H  . 

The  eigenvalues  k2r1  can  be  obtained  by   eliminating  Bn  from 
equation  (2.25) 


tan  k  H 

zn 


A    P. 


(2.53) 


Let 


X. 


k-H 

zn 


b 
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A 


(2.54) 


ton  Xn  =  -X„ 


bHVi^fJ 


(2.55) 


tan  Xn 

-   1 

Xn 

bHV^2 
c,2 

z 
CO 

H2 

(2.56) 


When  Cz  >  Cn  which  is  the  case  under  consideration, 
equation  (2.56)  possesses  real  roots  for  k„  and  Xn  ,  the 
later  varying  in  the  range  (n  -  1/2)  TT  1  X^£  n  <T/  .  On  the 
other  hand,  when  C2  <  C,  ,  the  roots  kn  are  complex  numbers 
with  negative  imaginary  parts.  In  the  later  case,  there- 
fore, the  factor  exp(-jkv;r)  in  equation  (2.52)  implies 
horizontal  attenuation.  In  the  case  of  a  fast  bottom,  there 
exist  solutions  which  suffer  no  horizontal  damping.  This 
condition  stems,  of  course,  from  the  fact  that  for  angles  of 
incidence  (measured  from  the  normal)  greater  than  critical 
angle,  no  power  is  transmitted  into  the  bottom.  Figures 
2.2,  2.3,  and  2.4  shows  the  z  dependence  of  the  first  three 
normal  modes  for  50  kHz,  120  kHz,  160  kHz  respectively. 
Figure  2.5  shows  the  r  dependence  of  the  sound  pressure. 
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Figure  2.2    Pressure  amplitude  vs  depth  for  first  mode 
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Figure  2.3    Pressure  amplitude  vs  depth  for  second  mode 
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Figure  2.4    Pressure  amplitude  vs  depth  for  third  mode. 
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C.   CUT  -  OFF  FREQUENCY 

Let  us  determine   the  cut-off  frequencies  of   the  normal 
modes.  It  follows  from  equation  (2.56)  that 


tonX„       _        z± 


\ 


bV^l^.X;        (2.57, 


It  is   clear  from  the  above   equation  that  the   nth  mode 
cannot  support  propagation  until 


"  (2.58) 


HU-i)  *v*-v*\ 


Therefore,   the  limiting  frequency  below  which  transmis- 
sion of  the  nth  mode  cannot  take  place  is  given  by 


f  (2n-1)c 

n       ~ 


4Hj~nr        (2-59> 


where    n   =    1,2,3, 

This  is  an  important  relation  between  the  layer  thick- 
ness H  and  the  cut-off  frequency  f  .  It  is  shown  in  Figure 
2.6  for  the  first  mode. 
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D.   PHASE  SPEED 

The  factor  exp  j  (cot-k^r+tf/4)  in  equation  (2.52)  allows 
an  obvious  interpretation  of  k^,  namely  kn  =  co/Cpn  where 
Cpn  is  the  phase  speed  of  the  nth  mode.  The  phase  speed  is 
the  speed  that  the  wave  front  (surface  of  constant  phase) 
propagates  in  r  direction.  The  phase  speed  at  the  cut-off 
frequency  has  the  value  of  the  speed  of  sound  in  the 
bottom  and  decreases  continuously  with  the  increasing 
frequency  towards  the  value  of  sound  speed  in  the  upper 
layer.   From  equation  (2.19)  we  can  write  [Ref.  3] 


2  2        2 

n  vcl)       zn 


OTJ   ~      zn  (2.60) 


2. 

Multiply   both    sides   by  H      gives 


2     2  2  2     2 

k  H      =r    (±M_\     _  k   H 

n  \     C,    /  zn  (2.61) 


or 


2.  2  2    2    2      2     2 

co  H      -     c,  (  k  H  +  k   H   ) 


zn 


or 


2    2  2*2 


(2.62) 
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where 


oL 


^    2 
k„H 


(2.63) 


Therefore  equation  (2.57)  can  be  written  as 


ton  X/7 
X„ 


-1/b 


HV^ 


w)  ~  (f ) 


(2.64) 


or 


tan   Xn 
X 


-   8 


VrtV-  X 


(2.65) 


where 


8    =      A±± 
A  ci 


(2.66) 


c/_  1 


Ct)- 


(2.67) 


From  equation  (2.65) 


*  =   XW1  +  ^  cot  (X  ) 
6  n 
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(2.68) 


We  know  that  C(n=      k^H   where  kn  =  co/C?n  . 
Cpr  is  the  phase  speed  of  the  nth  mode. 
Therefore 


C    = 


W 


•H 


oc. 


(2.69) 


From  equation  (2.62)  we  have 


coH   - 


2. 

i 

<* 

+ 

X 

n 

n 

(2.70) 


Hence 


=   C 


\ocnJ 


(2.71) 


Substituting  the  value  of  Xo/0(n  from  equation  (2.68)  we  get 


Pn 


R  ■■      **    1 

2 

1  +  6acot2(Xn). 

(2.72) 


(n-l/2)1Y  <  X„  <  n  1Y  for  the  nth  mode,  and  in  this  case 
k„  and  phase  speed  are  real  and  propagation  of  the  modes 
proceed  without  damping.  XM  is  a  function  of  angular 
frequency  and  the  frequency  corresponding  to  any  given 
value  of  Xn  can  be  obtained. 


33 


7     E       cj-H 
27TC, 


J± 


27TC 


2  7T£ 


f-62+  S^cotY 


(2.73) 


When      Xn  is    near    its   upper    limit    of   nff,then   cot   X    — ►  oO. 
Therefore    ft    is    very    large    and    from   equation    (2.72) 


Pn 


Hence  in  the  limit  of   very  high  frequency,   phase  speed 
approaches  the   sound  speed  CH . 
limits  ( n- 1/2 ) Tf then  cot  X  -»  0 


When  Xn  is  near   its  lower 


and 


C  =  C 

■pn      1 


i+e*]1  =  c, 


The  phase  speed,  therefore,  varies  from  the  value  Cz  at 
the  critical  frequency  down  to  C1  in  the  limit  of  very  high 
frequency.  Figure  (2.7)  shows  the  phase  speed  as  a  function 
of  frequency  for  first  three  modes. 

E.   GROUP  SPEED 

In  a  medium  in  which  the  phase  speed  varies  with 
frequency  and  if  the  pressure  pulse  from  a  point  source  is 
not  a  single  frequency  but  covers  a  moderately  broad 
spectrum,  then  the  presence  of  the  particular  frequency 
component  of  the  spectrum  at  a  large  range  would  occur  at  a 
time  associated  with  the  group  speed  (speed  with  which  the 
energy  associated  with  a  particular  frequency  will  travel). 
At  frequencies  co  for  which  the  phase  is  stationary, 
mutual    interference     will   be    negligible   and     these 
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frequencies  will  be  dominant  at  the  prescribed  values 
of  t  and  r.  The  values  of  t  and  r  for  which  the  phase 
is  stationary  for  a  given  frequency,  therefore,  determine 
the  rate  of  propagation  of  this  frequency  in  the  mutually 
interfering  train  of  sinusoidal  waves,  [Ref.  3],  called 
group  speed  C,  . 


f(w,r,t)      -       cot  -  kn(«)r   +    n  (2  74) 


but 


Therefore 


dl     =    t  -  r  d_ko_     =     o 


(2.75) 


%n   =    I    =    d« 

t  d k„  (2.76) 


co      =      k  c 

n    pn 


C n     =      dco_     -     d(knCPrl)  =        C     ,    k   dCp, 

*  dk,  dkn  r'  +    ndkn       (2.77) 


From  equation  (2.63)   we  have  k^  =  0(„/H  .   Substituting  this 
value  in  equation  (2.77),  we  get 
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can     =     c    +    o(  d  cP„ 

"dc*n 


9" 


Pin    '  /J 


(2.78) 


From   equation    (2.71)    we    can   write 


d  C„, 

d<*„ 


2 


2    -i 


+ 


(1) 


L   *„J 


d  Xn  Xn 


^  d*n 


<*„2 


(X    dcP„ 


d  Xp  _    Xn 

V°tf+  x; 


■iXnV 


(2.79) 


From  equation  (2.68)  we  get 


llX/i 
d«n 


f.  y  1  +  s2cot2  xff 

0(X„) 


(2.80) 


And  then  substitution  of  equation  (2.80)  and  equation  (2.79) 
into  equation  (2.78)  yields 


C, 


_2n   = 


1   +     e 


0(XJ 


(2.81) 


where 


^(Xn) 


2       2 

1  +  8  cot  X. 


SX„.cos  X,/sm3X^  ] 


(2.82) 
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Figure  (2.7)  shows  the  group  speed  as  a  function  of 
frequency  for  first  three  modes.  It  may  be  noted  that  group 
speed  is  always  less  than  the  phase  speed.  The  general 
features  for  each  mode  are,  [Ref.  2]. 

1 . The  first  arrival  would  be  nearly  sinusoidal  waves 
of  frequency   f n ,  where   fM  denotes   the   limiting 

frequency  for  the  n  th  mode.  The  wave  would  arrive 
at  t  =  r/C2  travelling  with  C9y,  =  C2  .  This  portion 
of  the  signal  is  called  ground  wave.  As  the  time 
progresses,  the  frequency  and  amplitude  in  the  ground 
wave  would  increase. 

2 . At  time  t  =  r/C,  a  new  train  of  high  frequency 
waves,  due  to  the  right  hand  branch  of  the  group 
velocity  curve,  would  suddenly  be  superimposed  on 
the  ground  wave  travelling  with  speed  C»n  =  C,  . 
This  new  high  frequency  wave  is  designated  as  water 
wave.  The  frequency  in  the  water  wave  would 
decrease  as  time  progresses  while  in  the  ground 
wave  the  frequency  would  continue  to  increase. 

3. Still  later  f^  and  fw  would  approach  each  other 
at  t=r/Cg0  ,  where  Cq0  denote  the  minimum  group  speed 
they  would  coincide.  The  pressure  then  would 
consists  of  a  single  frequency  fA  and  this  portion 
of  the  signal  is  called  Airy  phase.  The  sequence  of 
these  events  is  illustrated  in  Figure  2.8  . 


38 


-p 
o 


> 


4-> 
C 
0) 


a 

< 

□ 


z 
z 

— 
- 


CO 
C 

w 
en 

a 
-a 
> 

•H 

<1) 

u 

OS 


00 
CN 

(1) 

•H 


39 


Ill .  SOUND  PROPAGATION  IN  A  MODELED  SHALLOW  WATER 

ENVIRONMENT 

A.  INTRODUCTION 

The  purpose  of  this  experiment  was  to  investigate  exper- 
imentally the  behaviour  of  sound  propagation  in  a  laboratory 
facility  consisting  of  a  layer  of  water  overlaying  a  layer 
of  sand  and  to  compare  the  results  with  the  normal  mode 
theory. 

B.  EQUIPMENT 

1.  The  tank  used  for  experimental  work  is  made  of  alumi- 
nium plate  one  centimeter  thick.  The  interior  dimensions  of 
the  tank  are  : 

Length       =       9.95   meters 

Width        =       0.94   meters 

Depth  =  0.24  meters 
The  tank  has  a  carriage  which  can  be  used  to  adjust 
horizontal  and  vertical  positions  of  receiver.  The  rails  on 
which  the  carriage  moves  were  leveled  to  ensure  that 
receiver  remained  at  a  constant  depth  as  it  moved  along  the 
tank.  The  leveling  arrangement  was  not  precise  and  it  was 
found  that  the  total  vertical  variation  of  the  receiver 
was  within  +  1  mm.  All  the  joints  of  the  tank  were  sealed 
with  silicon  sealer  and  the  interior  of  the  tank  was  coated 
with  three  layers  of  primer  paint  (zinc  chromate)  to  avoid 
corrosion.  Wooden  wedges  are  fixed  on  the  inner  sides  of  the 
tank  to  minimize  the  reflection  of  sound.  The  tank  was 
filled  with  water  saturated  sand  to  form  an  11  cm  thick 
bottom  layer  and  a  thin  layer  of  water  (2  to  3  cm)  was  added 
above  the  sand.    The  sand  was  degassed  by  circulating  water 
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deep  into  the  sand  with  a  centrifugal  pump  at  moderate  pres- 
sure. This  process  was  repeated  several  times.  The  top 
surface  of  the  sand  was  leveled  with  reference  to  the  water 
surface.  This  leveled  surface  of  sand  was  horizonally  fiat 
within  +  1  mm.  The  layer  of  water  was  very  shallow  and  the 
thickness  of  this  layer  changes  with  time  due  to  evaporation 
of  the  water.  To  avoid  this,  the  tank  was  equiped  with  an 
automatic  water  level  control  which  maintains  the  desired 
level  of  water. 

2.  The  frequencies  used  in  this  work  ranged  from  25  kHz 
to  100  kHz.  The  source  transducer  was  built  from  a  rectan- 
gular crystal  of  piezoelectric  material  2.98  cm  long  and 
1.12  cm  wide.  The  directivity  index  of  the  transducer  at  50 
kHz  is  5.75  dB.  The  receiver  was  made  from  a  4  mm  by  4  mm 
piezoceramic  cylinder. 

C.   MEASUREMENT  OF  SPEED  OF  SOUND  AND  DENSITY 

Figure  3.1  shows  the  equipment  and  the  circuit  used  to 
measure  the  speed  of  sound  in  water-saturated  sand  and  in 
water.  Two  LC-5  transducers  were  used  as  source  and 
receiver.  A  high  frequency  tone  burst  of  the  order  of  one 
MHz  was  used  and  the  speed  of  sound  was  calculated  by 
direct  measurement  of  time-of-f light  from  the  oscilloscope 
and  the  distance  between  the  source  and  the  receiver. 
Experimentally  determined  values  for  the  speed  of  sound  in 
the  sand  and  in  water  are  shown  in  Table  I  and  Table  II.  The 
standard  deviation  of  the  measurements  is  an  indication  of 
the  precision  of  the  measurement  (approximately  1  percent). 
The  difference  between  the  measured  speed  of  sound  in  water 
and  the  value  predicted  by  theory  at  that  temperature  is  an 
indication  of  the  accuracy  of  the  measurement  (approximately 


0.95  percent).  The  density  of  the  sand  was  measured  by 
weighing  1000  cm  of  water-saturated  sand  and  calculating  its 
mass  The  measured  values  are  shown  in  Table  I . 
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TABLE  I 

Speed  o 

f  sound 

in 

sand 

No.          Time 

Division 

Scale 
(  As/div) 

Distance 
(cm) 

C 
(m/s) 

1           8.1 

50 

67.8 

1674.07 

2           6.4 

50 

54.3 

1696.88 

3           5.5 

50 

45.3 

1647.27 

4           4.6 

50 

38.0 

1652. 17 

5           6.3 

20 

21.  1 

1674. 6 

6           8.6 

10 

14.2 

1651. 16 

Average  Speed 

= 

1666.03 

m/s 

Standard  Deviati 

on     = 

17.58 

Percent  Deviation      = 

1.06 

% 

Measurement  For  Density  of  E 

and 

Mass  of  the  Container 

= 

542    gm 

Mass  of  1000  cm3  of 

sand  + 

Container   = 

2650.3  gm 

Density   =   2108.3 

kg/m 
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TABLE  II 

o 

Speed  of  sound  in  water  (Temperature  =  17   C) 


No, 

1 
2 
3 
4 
5 
6 
7 
8 


Time 
Division 

Scale 
(  lk  s/div) 

50 

Distance 
(cm) 

C 
(m/s) 

3.8 

27.5 

1447.3 

4.5 

50 

33.0 

1466.67 

6.0 

50 

43.5 

1450.0 

7.05 

50 

51.5 

1461.0 

8.8 

50 

64.5 

1465.91 

5.0 

100 

72.2 

1444.0 

5.7 

100 

83.7 

1468.42 

6.4 

100 

93.7 

1464.06 

Average  Speed  of  Sound      =  1458.43  m/s 
Standard  Deviation          =      9.1 

Speed  of  Sound  from  Theory  =  1472.39  m/s 

Percent  Deviation  =      0.95  % 
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D.   VOLTAGE  RESPONSE  AND  SENSITIVITY  OF  TRANSDUCERS 

The  basic  characteristics  of  the  source  and  receiver- 
were  not  known,  therefore  the  voltage  response  of  the  trans- 
mitter and  the  sensitivity  of  the  receiver  were  measured  by 
the  reciprocity  method  (Figure  3.2).  The  left  hand  side  of 
the  Figure  indicates  the  transducers  being  used  as  transmit- 
ters and  the  right  hand  side  as  receivers.  Experimentally 
measured  data  are  shown  in  Table  III  and  calculated  values 
of  voltage  response  and  sensitivity  are  shown  in  Table  IV. 
The  voltage  response  of  the  transmitter  is  plotted  in  Figure 
3.3  and  sensitivity  of  the  receiver  in  Figure  3.4  . 
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Figure  3.2    Calibration  of  transducers  by  reciprocity  method 
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TABLE  III 

Exp 
ca 

erimental  i 
librations 

measurements  for 
of  transducers 

Fre . 
(kHz) 

V* 
(volt) 

Vx 
(volt) 

la 

(amp) 

Vx 

(volt) 

20 

0.00046 

0.0015 

0.0051 

0.00036 

30 

0.00073 

0.002 

0.0075 

0.006 

40 

0.0012 

0.0016 

0.01 

0.012 

50 

0.00175 

0.0074 

0.0127 

0.034 

60 

0.00255 

0.0071 

0.0155 

0.044 

70 

0.003 

0.0227 

0.0185 

0.099 

80 

0.0046 

0.018 

0.0215 

0.224 

90 

0.006 

0.059 

0.025 

0.46 

100 

0.0075 

0.0625 

0.0275 

0.62 

110 

0.01 

0.185 

0.0305 

1.  15 

120 

0.0123 

0.089 

0.0315 

0.17 

130 

0.016 

0.085 

0.034 

0.193 

140 

0.015 

0.06 

0.037 

0.  1 

150 

0.0145 

0.093 

0.043 

0.21 

Distance 

between  source 

and  receiver   =   1 

m. 

Applied  volta 

ge 

=   5 

volts 
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TABLE 

IV 

Voltage  response  a 

nd 

sensitivity 

Fre 
(kHz 

Sensitivity 
)        dB(re  V/Pa) 

Voltage  Response 
dB(re  Pa/V  §  1  m) 

20 

-76.37 

5.91 

30 

-68.34 

0.38 

40 

-70.96 

1.06 

50 

-63.43 

6.84 

60 

-65.78 

8.83 

70 

-59.36 

12.50 

80 

-59.91 

11.03 

90 

-53.95 

15.39 

100 

-54.24 

16.18 

110 

-48.96 

20.32 

120 

-61.86 

26.87 

130 

-63.33 

27.94 

140 

-68.10 

29.69 

150 

-63.68 

29.07 
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E.   VALIDATION 

To  insure  that  the  sound  fields  set  up  in  the  tank  would 
be  good  approximations  of  the  sound  fields  in  the  absence  of 
sides  and  bottom  (e.g,  a  shallow  water  environment  of  infi- 
nite horizontal  extent),  measurement  were  made  on  the  sound 
field  set  up  in  the  water  and  in  the  sand  for  a  source  in 
the  water  layer. 

1 .   Sound  Pressure  as  a  Function  of  Receiver  Depth 

Figure  3.5  shows  the  circuit  arrangement  for  the 
measurement  of  the  variation  of  pressure  amplitude  with 
range  and  depth.  Continuous  waves  of  the  desired  frequency 
were  generated  by  a  frequency  generator  and  were  applied  to 
the  source  through  an  amplifier.  The  source  was  placed  at 
mid  depth  in  the  water  layer  which  was  2.7  cm  thick.  The 
signal  was  received  by  the  receiver  and  was  passed  through 
the  bandpass  filter  to  eliminate  noise.  The  output  from  the 
bandpass  filter  was  amplified  and  displayed  on  the  oscillo- 
scope. To  determine  the  variations  of  pressure  amplitude  as 
a  function  of  receiver  depth,  the  transmitter  was  kept  at  a 
constant  depth  in  the  water  layer  while  the  measurements 
were  made  with  the  receiver  at  a  fixed  range  beginning  from 
the  top  surface  of  the  water  layer  down  into  the  bottom 
layer  in  steps  of  0.5  centimeter.  Three  frequencies  between 
the  cut-offs  of  two  first  modes  were  selected  for  measure- 
ments. The  data  are  tabulated  in  Tables  V,  VI,  and  VII  and 
plotted  in  Figures  3.6,  3.7,  and  3.8  .  The  experimental  and 
theoretical  results  agree  satisfactorily  within  the  water 
layer.  Near  the  cut-off  frequency  of  the  propagating  mode, 
the  experimental  amplitude  appears  to  decay  more  rapidly 
with  depth  in  the  bottom  layer  than  theory  predict.  Decay  is 
slower  for  intermediate  and  high  frequencies. 
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TABLE 

V 

Sound 
f 

pressure 
=35  kHz, 

as 

n 

function 
=  1,  L   = 

of  depth: 
1.5  m 

No. 

Depth ( cm) 

Volts 

Normalized 

1 

0.0 

0.012 

0.12 

2 

0.5 

0.028 

0.28 

3 

1.0 

0.056 

0.56 

4 

1.5 

0.081 

0.81 

5 

2.0 

0.097 

0.97 

6 

2.5 

0.1 

1.0 

7 

3.0 

0.097 

0.97 

8 

3.5 

0.083 

0.83 

9 

4.0 

0.067 

0.67 

10 

4.5 

0.052 

0.52 

11 

5.0 

0.042 

0.42 

12 

5.5 

0.041 

0.41 

13 

6.0 

0.038 

0.38 

14 

6.5 

0.034 

0.34 

15 

7.0 

0.028 

0.28 

16 

7.5 

0.023 

0.23 

17 

8.0 

0.018 

0.  18 

18 

8.5 

0.015 

0.  15 

19 

9.0 

0.016 

0.  16 
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Pressure  variation  as  function  of  receiver 
depth: f  =  35  kHz  n  =  1. 
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TABLE 

VI 

Sound 
f 

P 

ressure 
50  kHz, 

as 

n 

function 
=  1,  L  = 

of  depth: 
1.5  m 

No. 

D 

spth( cm) 

Volts 

Normalized 

1 

0.0 

0.04 

0.06 

2 

0.5 

0.21 

0.31 

3 

1.0 

0.41 

0.61 

4 

1.5 

0.58 

0.87 

5 

2.0 

0.67 

1.0 

6 

2.5 

0.65 

0.97 

7 

3.0 

0.58 

0.87 

8 

3.5 

0.435 

0.65 

9 

4.0 

0.34 

0.51 

10 

4.5 

0.245 

0.37 

11 

5.0 

0.  15 

0.22 

12 

5.5 

0.1 

0.  15 

13 

6.0 

0.063 

0.09 

14 

6.5 

0.05 

0.07 

15 

7.0 

0.05 

0.07 

16 

7.5 

0.047 

0.07 

17 

8.0 

0.027 

0.04 

18 

8.5 

0.02 

0.03 

19 

9.0 

0.04 

0.06 
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Pressure  variation  as  function  of  receiver 
depth :f  =  50  kHz  n  =  1. 
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TABLE 

VII 

Sound  pressure 
f  =  70  kHz, 

as 

n 

function 
=  1,  L  = 

of  depth: 
2.0  m 

No. 

Depth ( 

cm) 

Volts 

Normalized 

1 

0.0 

0.095 

0.05 

2 

0.5 

0.53 

0.28 

3 

1.0 

1.15 

0.6 

4 

1.5 

1.7 

0.89 

5 

2.0 

1.9 

1.0 

6 

2.5 

1.75 

0.92 

7 

3.0 

1.2 

0.63 

8 

3.5 

0.7 

0.37 

9 

4.0 

0.24 

0.13 

10 

4.5 

0.205 

0.11 

11 

5.0 

0.2 

0.  1 

12 

5.5 

0.  195 

°  0.1 

13 

6.0 

0.  13 

0.07 

14 

6.5 

0.043 

0.02 

15 

7.0 

0.06 

0.03 

16 

7.5 

0.065 

0.03 

17 

8.0 

0.1 

0.05 

18 

8.5 

0.  105 

0.055 

19 

9.0 

0.073 

0.04 
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Pressure  variation  as  function  of  receiver 
depth: f  =  70  kHz  n  =  1. 
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2 .   Sound  Pressure  as  a  Function  of  Range 

To  measure  the  variation  of  peak  pressure  amplitude 
with  range,  the  transmitter  was  kept  at  constant  depth  and 
the  receiver's  depth  was  adjusted  slightly  at  each  range  to 
measure  peak  amplitude  at  that  range.  Precautions  had  to  be 
taken  to  ensure  that  the  water  surface  is  disturbed  as 
little  as  possible  because  variation  in  water  depth  had  a 
significant  effect  on  the  pressure  amplitude.  Measurements 
were  made  for  three  frequencies  between  the  cut-offs  of 
the  first  two  modes.  The  data  are  tabulated  in  Tables  VIII, 
IX,  and  X  and  is  plotted  in  Figures  3.9,  3.10,  and  3.11 
respectively.  Assuming  a  simple  model  for  which  the  loss  in 
the  bottom  effects  that  portion  of  the  energy  density  of  the 
normal  mode  found  in  the  bottom,  which  then  "robs"  the 
energy  density  in  the  water  layer,  then  the  effective 
absorption  coefficient  can  be  estimated  as  (Ref.  personal 
communication  with  Prof.  Alan  B.  Coppens) 

a  *  ±  F  Aj>      dB/m 
4   A         7 

Where  Ab/A  is  fraction  of  the  energy  density  of  the  normal 
mode  in  the  bottom  and  1/4F  (F  is  the  frequency  in  kHz) 
yields  the  expected  attenuation  to  be  found  in  the  sand  used 
for  the  bottom  [Ref.  4].  Comparison  of  the  measured  data 
with  the  simplestic  model  is  very  good,  suggesting  that 
indeed  the  sound  field  is  subjected  to  minimal  effects  from 
the  side  walls. 
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TABLE 

;  viii 

Sound 

HP 

ressure  as 
=  2  cm  f  = 

funct 
50  kH 

ion  of 
z  n  =  1 

range 

No. 

D 

istance(m) 

Volts 

No 

rmalized 

1 

1.00 

0.45 

1.0 

2 

1.25 

0.34 

0.756 

3 

1.50 

0.255 

0.567 

4 

1.75 

0.  195 

0.433 

5 

2.00 

0.2 

0.444 

6 

2.25 

0.  168 

0.373 

7 

2.50 

0.  108 

0.24 

8 

2.75 

0.  155 

0.344 

9 

3.00 

0.  105 

0.233 

10 

3.25 

0.078 

0.  173 

11 

3.50 

0.022 

0.049 

12 

3.75 

0.047 

0.104 

13 

4.00 

0.059 

0.  131 

14 

4.25 

0.057 

0.  127 

15 

4.50 

0.043 

0.095 

16 

4.75 

0.042 

0.093 

17 

5.00 

0.031 

0.069 
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TABLE 

IX 

Sound 

pressure 
H  =  2  cm 

as 

f  = 

fu 
70 

net 
kH 

ion  of 
z  n  =  1 

range 

; 

No. 

Distance(m) 

Volts 

No 

rmalized 

1 

1.00 

2.08 

1.0 

2 

1.25 

1.6 

0.769 

3 

1.50 

1.7 

0.817 

4 

1.75 

1.0 

0.481 

5 

2.00 

0.93 

0.447 

6 

2.25 

1.1 

0.529 

7 

2.50 

0.85 

0.409 

8 

2.75 

0.68 

0.327 

9 

3.00 

0.87 

0.418 

10 

3.25 

0.62 

0.298 

11 

3.50 

0.17 

0.082 

12 
13 
14 
15 
16 
17 

3.75 
4.00 
4.25 
4.50 
4.75 
5.00 

0.61 

0.58 

0.36 

0.315 

0.3 

0.29 

0.293 
0.279 
0.173 
0.  151 
0.  144 
0.  139 
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TABLE 

X 

Sound 

HP 

ressure  as 

=  2  cm  f  = 

func 
100 

:tion  of 

kHz  n  = 

range 

: 

No. 

D 

istance(m) 

Volts 

No 

rmalized 

1 

1.00 

1.95 

1.0 

2 

1.25 

1.4 

0.718 

3 

1.50 

1.4 

0.718 

4 

1.75 

1.05 

0.538 

5 

2.00 

1.0 

0.513 

6 

2.25 

0.88 

0.451 

7 

2.50 

0.83 

0.42  6 

8 

2.75 

0.85 

0.436 

9 

3.00 

0.74 

0.379 

10 

3.25 

0.63 

0.323 

11 

3.50 

0.7 

0.359 

12 

3.75 

0.28 

0.  144 

13 

4.00 

0.65 

0.333 

14 

4.25 

0.5 

0.256 

15 

4.50 

0.41 

0.21 

16 

4.75 

0.49 

0.251 

17 

5.00 

0.41 

0.21 
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3 .   Sound   Pressure  Variation   Across  the   Width  of   the 
Tank 

Sound  pressure  variations  were  measured  across  the 
width  of  the  tank  at  a  one  meter  distance  from  the  source 
for  three  different  frequencies.  The  measured  data  shown  in 
Table  XI,  XII,  and  XIII  and  plotted  in  Figures  3.12,  3.13, 
and  3.14.  The  results  show  slight  effects  of  interference 
either  from  the  walls  or  from  the  bottom  of  the  sand  layer. 
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TABLE 

XI 

Transverse 

measurement 
f  =  50  kHz, 

of 
L 

sound  pi 

=  1  m 

ressure : 

Distance 
(cm) 

Volts 

Distance 
(cm) 

Volts 

-32 

0.38 

-3.5 

0.4 

-30.5 

0.215 

-1.7 

0.55 

-29 

0.38 

0 

0.40 

-27.5 

0.34 

2 

0.55 

-26.4 

0.39 

4 

0.4 

-25.2 

0.32 

6.7 

0.53 

-23.7 

0.47 

8.5 

0.34 

-22 

0.33 

10.7 

0.49 

-20 

0.5 

13.8 

0.36 

-18 

0.33 

16.3 

0.49 

-16.3 

0.5 

18.8 

0.32 

-15 

0.42 

20.2 

0.4 

-13.5 

0.5 

22.2 

0.31 

-11.8 

0.36 

23.5 

0.38 

-9.9 

0.53 

24.7 

0.33 

-7.2 

0.42 

26 

0.37 

-5.3 

0.53 

27.7 

0.25 
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TABLE 

XII 

Transverse 

measurement 
f  =  70  kHz, 

of  sound  p: 
L  =  1  m 

ressure : 

Distance 
(cm) 

Volts 

Distance 
(cm) 

Volts 

-32 

1.7 

-1.3 

1.7 

-30.6 

1.25 

0.4 

2.47 

-29.4 

1.4 

2 

1.6 

-28.6 

1.25 

3.7 

2.  17 

-27.2 

1.8 

5.2 

1.35 

-26 

1.25 

7 

2.23 

-24.9 

1.75 

8.7 

1.45 

-23.7 

1.2 

10 

2.  12 

-22.5 

2.05 

11 

1.85 

-21.3 

1.2 

12.3 

2.2 

-19.8 

2.2 

13.7 

1.4 

-18.4 

1.5 

15 

2.  1 

-16.8 

1.9 

16.5 

1.65 

-15.  6 

1.7 

17.8 

1.9 

-14.2 

2.1 

18.8 

1.57 

-13.3 

1.82 

19.8 

1.75 

-12 

1.98 

21.2 

1.3 

-10.7 

1.75 

22.5 

1.  62 

-9.2 

2.2 

23.6 

1.35 

-7.8 

1.55 

24.5 

1.6 

-6 

2.5 

26 

1.0 

-4.5 

1.65 

27.8 

1.5 

-3 

2.32 
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TABLE 

XIII 

Transverse 

measurement 
f=  100  kHz, 

of  sound  p: 
L  =  1  m 

ressure : 

Distance 
(cm) 

Volts 

Distance 
(cm) 

Volts 

-32 

1.25 

-4 

1.47 

-30.8 

0.82 

-1.9 

1.9 

-29.6 

1.5 

-0.4 

1.52 

-28.4 

0.94 

0.5 

1.7 

-27 

1.4 

1.5 

1.6 

-25.5 

0.94 

2.4 

1.7 

-23  .8 

1.62 

3.4 

1.62 

-22.6 

1.15 

5.3 

1.22 

-20.5 

1.76 

7.8 

1.67 

-18.7 

1.55 

10 

1.  1 

-17.5 

1.8 

12.5 

1.35 

-16.4 

1.5 

15.6 

0.82 

-15.4 

1.9 

16.5 

1.22 

-14.4 

1.5 

18 

0.63 

-13.5 

1.8 

19.8 

1.05 

-13 

1.6 

21.2 

0.67 

-11.7 

2.2 

22.5 

0.95 

-10.7 

1.65 

24.7 

0.49 

-9.5 

1.8 

26 

0.68 

-8.7 

1.65 

27.5 

0.42 

-7 

1.88 
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F.   SUMMARY  AND  CONCLUSION 

In  the  preceding  sections,  the  experimental  curves  for 
the  pressure  amplitude  as  a  function  of  depth  and  range  for 
three  different  frequencies  for  which  only  the  first  normal 
mode  propagates  have  been  compared  with  normal  mode  theory 
and  a  simple  model  for  attenuation  of  the  normal  mode.  The 
depth  dependent  curves  both  in  the  water  and  in  the  sand 
agrees  satisfactorily  with  the  theory.  The  range  dependence 
obtained  in  the  water  exhibit  excellent  agreement  with  the 
simple  modeling  of  losses  resulting  from  the  relatively 
highly  attenuating  bottom.  Measurements  both  in  range  and 
transverse  dimensions  suggests  that  the  sound  field  is 
subjected  to  minimal  effects  from  the  side  walls.  The 
leveling  arrangement  needs  improvement  because  at  higher 
frequencies  minor  alterations  in  water  depth  changes  the 
pressure  amplitude  considerably. 
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